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(D ■ Abstract 

'""'■ We study gravitational and electromagnetic perturbation around the squashed Kaluza-Klein 

^ . black holes with charge. Since the black hole spacetime focused on this paper have SU{2) x U{1) ~ 

\^ \ U{2) symmetry, we can separate the variables of the equations for perturbations by using Wigner 

T-H ' function -D^Af which is the irreducible representation of the symmetry. In this paper, we mainly 

f— ^ ■ treat J = modes which preserve SU{2) symmetry. We derive the master equations for the J = 

o : 

,—1 , modes and discuss the stability of these modes. We show that the modes of J = and K = Q, ±2 

.5^ ' and the modes oi K = ib(J + 2) are stable against small perturbations from the positivity of the 

;^ ■ effective potential. As for J = 0,i^ = ±1 modes, since there are domains where the effective 

potential is negative except for maximally charged case, it is hard to show the stability of these 
modes in general. To show stability for J = 0,i^ = ±1 modes in general is open issue. However, 
we can show the stability for J = 0, X = ±1 modes in maximally charged case where the effective 
potential are positive out side of the horizon. 
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I. INTRODUCTION 



In recent years, the studies on higher dimensional black holes have attracted much atten- 
tion in the context of the brane world scenario. One of the main reason is that mini black 
hole creation in a particle collider was suggested [iHsl based on some brane world model. 
Such physical phenomena are expected to give us a piece of evidence for the existence of 
extra dimensions and to draw some information towards quantum gravity. 

It is known that higher-dimensional black objects admit various topology of horizon unlike 

ligher 



16|. 



four- dimensional case. In fact, there are many exotic solution of Einstein equation in 
dimensions in addition to black holes, e.g. black rings, strings, Starns and di-rings [7 
To predict what type of black objects are created in a particle collider, one of the important 
criterion is their stability against small gravitational perturbations. In |l7|, Kodama and 
Ishibashi showed that the equation of motion for perturbation of the maximally symmetric 
higher-dimensional black holes reduces to single master equation for each mode. They also 
showed the stability of the higher- dimensional Schwarzschild black holes. Recently, much 



effort has been devoted to reveal the stabi 
complicated structure by many authors 118 



ity of higher-dimensional black holes with more 



26|. 



From a realistic point of view, the extra dimensions need to be compactified to recon- 
cile the higher- dimensional gravity theory with our apparently four-dimensional world. We 
call the higher-dimensional black holes on the spacetime with compact extra dimensions 
Kaluza-Klein black holes. In general, it is difficult to construct exact solutions describing 
Kaluza-Klein black holes because of the less symmetry than the asymptotically flat case. 
However, if we consider the spacetime with twisted extra-dimensions, we can construct 



28 



in the class 



such exact solutions, i.e. squashed Kaluza-Klein (SqKK) black holes [27, 
of cohomogeneity-one symmetry. The topology of the horizon of this SqKK black holes is 
5*^, while it looks like four- dimensional black holes with a circle as an internal space in the 
asymptotic region. 

Recently, much effort has been devoted to reveal the properties of squashed Kaluza- 



Klein black holes. In 



29I47 



47| generalizations of SqKK black holes are studied. Several 



aspects of S qKK black holes are also discussed, e.g. thermodynamics |48l-l50|. Hawking 
radiation 5ll-l53|. gravitational collapse 5J], behavior of Klein-Gordon equation [55], stabil- 



ity for neutral case |56|^, quasinormal modes [57|, |58|, geodetic precession [59|, Kerr/CFT 
correspondence 60| and gravitational lensing effects 61 |. 

In this paper, we study gravitational and electromagnetic perturbation around the 
squashed Kaluza-Klein black holes with charge 3l| by extending the analysis in [56] . SqKK 



black hole with charge have t 
the same technics used in 120 



le SU{2) X f/(l) symmetry as in the neutral case, so we can use 



56 



62| to analyze the perturbations of SqKK black holes with 



charge. We expand the perturbation variables in terms of Wigner function -D^Af ? which is 
the irreducible representation of the symmetry characterized by three indices J, M, K. Since 
the instability empirically appears in the lower modes, we mainly focus on the perturbations 
in J = modes which preserve SU{2) symmetry. In this paper, we derive the master equa- 
tions for {J = 0,M = 0,K = 0, ±1, ±2) modes which have SU{2) symmetry and the highest 
modes {K = ±(J + 2)). Using the effective potential functions in the master equations, we 
discuss the stability for these modes. 

There are several motivations to study the perturbation for SqKK black holes with charge. 
One is to understand the property of the spacetime deeper by using small perturbation as 



a probe such as quasinormal modes. Next, the stability of Sc 
the arguments in the physics around the black holes 48l-l53 



KK b 



57|-l5S, 



ack holes is needed for 



61| to be meaningful. 



Finally, it is also useful to consider the perturbation for SqKK black holes with charge to 
understand the effective theories reduced from unified theories based on string theory which 
usually contains not only gravity but also gauge fields. 

The organization of this paper is as follows. In section [TTl we review the geometry of 
SqKK black holes with charge and the formalism to classify metric perturbations based 
on the symmetry. In section IIIIl we derive the master equations for master variables. By 
analyzing these equations, we discuss the stability of SqKK black holes with charge. The 
final section is devoted to the discussion. 



^ Note that the stabihty of J = 0,K = ±1 modes were not shown analytically since there are some 
mistakes in the discussion in J = 0,iir = ±1 modes in 56|. It is hard to find S-deformation function 
to show stability for these modes analytically. To find S-deformation function analytically is open issue. 
However, the stability for these modes have been confirmed by numerical method in 57 1. 



II. SQUASHED KALUZA-KLEIN BLACK HOLES WITH CHARGE 



In this section, we review the geometry of SqKK black holes with charge 31| and the 
formalism to classify metric perturbations based on the symmetry 20|, |56|, l62 |. 



A. Geometry of squashed Kaluza-Klein black holes with charge 

We start with the five-dimensional Einstein-Maxwell system whose action is given by 

^=1^1 ^'"^^^ ^^ " ^"''^'^^ ' ^^^ 

where G, R and F^j^^, = df^A^ — d^A^ are the five- dimensional gravitational constant, the 
Ricci scalar curvature and the Maxwell field strength with the gauge potential A^. From 
this action, we obtain the Einstein equation 



^flU n 9fJ,U '^-'■fiui 



(2) 



with 



and the Maxwell equation 



-'- ^u ^ pbX^ u .Qfiu^ aj3^ 



V^F^^ = 0. 



a/3 



(3) 



(4) 



SqKK black holes with charge [3l| are solutions of Eqs. 02]) and ^ whose metric and 
gauge potential are given by 

ds' = -F{p)de + ^dp' + p'K{pf[{ay + [a^f] + ^^° + ^^^)^^^; + ^'^ (a^)^ (5) 



np) 

^ 2 p ' 

where the functions F{p) and K{p) are 



K{py 



(6) 



np) 



{p - p+){p - p') 
p" 



P 



(7) 



and po and p± are constants which satisfy 



P+> P-> 0, p- + Po > 0, 



(8) 



and the invariant one forms a"" (a = 1, 2, 3) of SU{2) are given by 

a^ = — sin if) dO + cos i}) sin Odcj)^ a'^ = cos ipdO + sin ip sin 9 d(f), a^ = dip + cos 9 dcp. (9) 

The domains of angular coordinates are O<^<7r,O<0< 2tt, < ip < Air, and radial 
coordinate p runs in the range < p < oo. Horizons locate at p = p±. 
In the region far from the horizon, the metric (j5]) becomes 

2 

ds'^ ~ -dt^ + dp^ + p^{d9'^ + sin^ 9d(t)'^) + ^(# + cos9d(t)f + C(l/p), (10) 



where r^o = 2^J{pJ^ + po){p- + po) is the size of the extra-dimension. So we can see that the 
spacetime behaves effectively four-dimensional Minkowski spacetime if we focus on physical 
phenomena whose typical size is much larger than r^. On the other hand, in the vicinity 
of the horizon, the metric ([5]) has five-dimensional nature. To see this we introduce new 
coordinates rj, r and new parameters r± as 

, 4p2 por2 Pari 

t = ^-V, P=- 2^ P± = -2 2"' (11) 



(y ^ (V" ^ (y ^ (V» ■'■ 



then the metric (151) becomes 



ds^ = -f(r)dv' + ^dr^ + r' [k{r) [{a^f + [a^f] + [a'f] , (12) 



with 



/w^ "-"'-f "-' . >=ix)- "--;^^^-'-' . (13) 

In this coordinate, the horizons locate at r = rj. and the shape of horizon is a squashed 5"^ 
which is characterized by the function k{r). Especially, if the size of the horizon is much 
smaller than that of extra-dimension Too ^ ^±, then we can see that the function k{r) ~ 1 
and the metric f ll2p behaves almost five-dimensional Reissner-Nordstrom black holes in the 
vicinity of the horizon. 



The Komar mass M and electric charge Q for this SqKK black hole are given by 3l|, |49| 



M=^-^{p^ + P^l Q = ^^Vp^. (14) 

For later calculations, it is convenient to define the new invariant forms, 

a^ = l{a'T^cr'). (15) 

The metric ([5]) can be rewritten by using a"^, 

ds' = -Fipidt' + §M!v + 4p'K{pfa-a- + (^*±M«±ih) (,3)^ . (le) 



B. Classification of the perturbations based on the symmetry 



Next, we introduce the formalism to classify metric perturbations based on the symme- 
try 20|, l56|, |62|]. The metric i^ has spatial symmetry of SU{2) x U{1) c^ U{2) which is 
generated by the Killing vectors 



^x 



COS( 



H — : — -Onh — cot u sm 



sin^ 
cos I 



4 = - sin (j)de + ^—p;^i, 
sm fc' 

63 = d^. 



d,i, — cot 9 cos 



(17) 



We treat the metric perturbation g^^, + h^^ and the electromagnetic perturbation A^ + 5A^, 
where Qfj,,, and A^ are the background quantities. Because of the background symmetry, 
these perturbations can be expanded by using Wigner functions which are known as the 
irreducible representation of SU{2) x U{1) [20|, |56|, l62|. Scalar Wigner function is defined by 



L^Di^, = J{J+l)Di,,, L-^Di^ = MDi^, W,Di^ = KD 



KM ? 



(18) 



where La := i^a (a = x, y, z), W3 := ie^, L^ := L^ + Ly + L^, and J, K, M are integers or 
half integers^ satisfying J > 0, \K\ < J, \M\ < J. Vector Wigner functions are constructed 
by the action of invariant forms a"" on scalar Wigner functions. 



DIk 



atDK-i, i\K-l\<J), 



crtDK, 



{\K\<J). 



(19) 



Here, we have omitted the subscript J, M. The vector Wigner functions satisfy 



L'Dl^ = JiJ + l)Dl^, L,DIj, = MDIj,, W,Dl^ = KDl^, 



(20) 



K, M take integers iff J takes integers. 



where a = ±, 3 and operations are defined by Lie derivatives; W^D\ p^ := CwzD\j^ and 
LaDfj^ := Clo^D^k- Tensor Wigner functions D^^j^ are defined by 



D, 



++ 

ij,K 



^M^^-2 (|i^-2|<J), 



DSk 
r)33 



ataTDK {\K\ < J), 



atapK-i {\K-l\<J), 

arajDK+2 {\K + 2\<J), 

a7apK+i {\K + 1\<J), 

afapK i\K\ < J). 



(21) 



One can check that D"^^^ forms the irreducible representation: 



■\ab 



L'DtlK = AJ + l)DS,K, ^.I^g;, = MDg^, W,D^,^^K 






(22) 



The tensor field h^i^ can be divided into three parts, Hab, hAi, hij, {A, B = t, p) which behave 
as scalar, vector and tensor within the submanifold 6,(f),ilj. Similarly, the vector field 6Afj, 
can be divided into two parts, 6Aa and 6Ai. These perturbations can be expanded using 
the Winger functions as 



1 1 iii/LttX' LtiAj 



dx^dx^ J2 h%{x^)DK{x') + 2dx^dx' Y^ hf^{x'^)Dlj^{x' 

K K 



-rdx^dx^y Xbi^WA^ 



K 



5Ai,dx^ 



dx^Y5A'i{x^)DK{x' 

K 



dx'Y5Aa{x'^)Dlj^{x' 

K 



(23) 



Because perturbed quantities are expanded in terms of the representation of the spacetime 
symmetry, no coupling appears between coefficients with different sets of indices (J, M, K) 
in the perturbed equations. In addition, we utilize the Fourier expansion with respect to 
the time coordinate t. 

Interestingly, without explicit calculation, we can reveal the structure of couplings be- 
tween coefficients with the same (J, M, K). First, since the index K is shifted in the defini- 
tion of vector and tensor harmonics, then the coefficients h^^, h^^ and /i^ exist only if K 
satisfies the inequality listed in the following table: 



K+ 


hA+,h+3,5A+ 


hAB, hA3, h^-, /i33, 5Aa, M3 


hA^,h_s,SA_ 


h__ 


K-2 < J 


\K-1\ < J 


K\<J 


\K+1 < J 


\K + 2\ < J 



Therefore, for zero mode J = 0, we can classify the coefficients by possible K as follows: 
J = 0; 



/^++ 


hA+,h+3,SA+ 


fiAB, hA3, h+-, /i33, 6 A A, 6 As 


hA-,h-3,6A_ 


/i__ 


K = 2 












K = 1 












K = 












K = -l 












K = -2 



Apparently, for hj^^ and h , we can obtain equations for a single variable, respectively. 

For other sets of coefficients {hA+, h^^, 5^+), (^ab? ^A3) h+-, ^33) ^Aa, SA3), {Ha-, /i~3, SA^) 
they are coupled to the coefficients in each set. As we will see later, after fixing the gauge 
symmetry, we have the coupled master equations for one or two master variables in each set. 

We can also easily see that /i++ in (J, M, K = J+2) modes and h in (J, M, K = — (J + 2)) 

modes are always decoupled. The perturbed equations for these modes can be reduced to 
the master equations for the single variables, respectively. We discuss the stability for J = 
modes and K = ±( J + 2) modes in the next section. 

III. STABILITY ANALYSIS 

The linearized Einstein equation is 



SGau - '25Ta^ = 0, 



(24) 



where the linearized Einstein tensor 6G^i, and the linearized stress tensor 6Tnu are given by 



where 



6G 



HV 



ST, 



^u 



-^M.(V^V^/ipA - Wph + h'^'R.x) - K.R], 
+6F,^F^^ - ^g^JF^^F'^", 



(25) 



(26) 



h = g'^'h^,, 5F^, = d,{5A,) - <9,(M^), 



(27) 



and R^u, F^^ are the background quantities. Note that V^ denotes the covariant derivative 
with respect to the background metric g^j_u- The Maxwell equation for the perturbations is 

= 0. (28) 

We treat two kind of gauge transformation. One is related to infinitesimal coordinate 
transformation : 

V -^ V + V^e^ + V^^M. (29) 

6A^ ^ 6A^ + rV.A^ + A.V^r , (30) 

where ^^ is arbitrary vector field. The other is related to f/(l) gauge: 

SA^ -^ SAf, + V^x, (31) 

where x is arbitrary scalar field. 

A. Zero mode perturbations ( J = 0, M = 0) 

In the case of J = 0,M = 0, there are five modes, K = ±2, ±1,0. We analyze these 
modes separately. These modes correspond to perturbations with SU{2) symmetry. 

1. K = ±2 modes 

hi K = ±2 modes, there are two coefficients /i+4. and h , and these are gauge invariant. 

Because the perturbed metric h^^{x'^)dx^dx'^ is real, h must be the complex conjugate of 

/i++, i.e. /i++ = h where the bar denotes the complex conjugate. So, it is sufficient to 

treat only /i++. We set /i^^ as 

h^^{x^)dx^'dx'' = /i++(p)e-^'^V+a+. (32) 

Inserting Eq. (|32l) into (++) component of Eq. ( !24l) . we obtain the perturbation equation 
for /i+_|. whose explicit form is written in Appendix |Al By introducing the master variable 



and switching to the tortoise coordinate p* defined by 

dp^ _ K{p) 
dp ~ F{p) ' 

we get the perturbation equation in the Schrodinger form 

rf2 



(34) 



2*^2 + V2{p)^2 = ^^$2, (35) 



where the potential V2(p) is defined by 

{p - P+){p - P+ + P- 



V2ip) 



16po(po + P-)p^(Po + P-+ p- p+)^ 

X [4p+(3p_ + 8p+)p^ + 4p_p+(9p_ + 41p+)pj^ + 4pV+(9p- + 74p+)p2 

+12pip+(p_ + 19p+)po + 64pip^ + ((88p+ - 9p_)p^ 

+ {-27 pi + 418p+p_ + 200pI)pI + P-{-27pl + 700p+p_ + 655p^)p^ 

+ (-9pi + 498p+p3_ + 711p^p^_)po + 128pip+(p_ + 2p+))(p - p+) 

+ (35p^ + (149p„ + 442p+)p3 + (257p^ + 1400p+p_ + 355p^)p2 

+p_(207p^ + 1470p+p_ + 739p^)po + 64p^(p^ + 8p+/5_ + 6p^))(p - p+f 

+ (200p|^ + 8(80p_ + 91p+)p^ + 8(87p^ + 187p+p_ + 32p^)po 
+256p4p2_ ^ 3^^^~_ ^ ^^))(^ _ ^^)3 

+ (352p2 + (736p_ + 512p+)po + 64(6p2_ + 8p+p_ + pl))ip - p+f 

+ (256po + 128(2p_ + p+))(p - p+f + 64(p - p^f] . (36) 

Here, we introduced new parameters p_ and po defined by 

P- := P+ - P-, Po := Po + P-- (37) 

From Eq. (jHj), these parameters satisfy the following inequalities, 

< p- < p+, Po > 0. (38) 

Since the terms in square bracket in Eq. (!36|) are given in a power series expansion of (p — p^) 
and the expansion coefficients of these terms are positive, we can see V2 > in the region 
p+ < p < 00. Typical profiles of the potential V2 are plotted in FigsHJISl One should be 
noted that we have equality y/SQ / {2M) = 1 in maximally charged case p+ = p_ from fin|) . 
and we normalized the depicted effective potential by using the size of extra dimension. The 
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FIG. 1: The effective potential r'^V2 in maximally charged case. 
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FIG. 2: The effective potential r^F2 in V3Q/2M = 2/3 case. 
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FIG. 3: The effective potential r^Fs in V3Q/2M = 1/3 case. 
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desired solution for $2 is square integrable in the domain —00 < p* < 00 with real uj^. 
Multiplying both sides of Eq. fl35|) by <l'2, and integrating it, we obtain 

2 



rfp* 



(i<l>2 



dp^ 



+ V"2|$2| 



7 ~\ 00 

rfp* 



u? I c?p* |<l>2|^ 



(39) 



Because the boundary term vanishes, the positivity of V2 means w^ > 0. Therefore, We 
conclude that the background spacetime is stable with respect to the K = ±2 perturbations. 

2. K = ±1 modes 



Because of the relations hA+ = h-A-, Z^+s = /i-3 and 6A^ = 6A_, it is sufficient to consider 
only hA+, /i3+ and 5^4+ . Thus, we assume 



h^^dx^dx" = 2hA+{p)e'''''dx^a+ + 2/i+3(p)e-*"V+(T 



(40) 



SA^dx" = 5AMe (T 



(41) 



Here, we set the gauge ^^dx^ as 



i,dx^ = i+{p)e''-'a 



(42) 



and under the gauge transformations ( !29l) . the metric perturbations transform as 



ht+ -^ ht+ - iuj^+, 

2p + Po 



hp+ — )■ /ip+ 



/ls_L -)■ /l.q_L - 



e+ + 



d^, 



Pip + Po) ^ ' dp 

i (p2 + 2pop - P-P+ - (p- 



P+)Po 



*3+ 



■3+ 



So we can impose the gauge condition 



^3+ 



(p + poy 



0, 



(43) 
(44) 

(45) 



(46) 



which completely fixes the gauge freedom. Substituting Eqs. fj40|) . fj4T]) and fj46|) into Eqs. 
f l2^ and (!28|) . we obtain perturbation equations whose explicit forms are given in Appendix 
lAl Eliminating ht+ from these equations, and defining new variables 

, (P - P+){p - P^){p+P- + Po(p+ + P- - 2p) - p^) 



>IG 



p"'{p + Pof/' 



-h 



P+' 



(47) 
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^,, := -^/-^. (48) 

we have the coupled master equations for the gravitational and the electromagnetic pertur- 
bations 



- ~pi^^G + Vn{p)(j)iG + Vi2{p)(j)lE = W 01G, (49) 

- 'T^^^E + V22{p)(plE + V2i{p)(piG = Oj'^(plE, (50) 



where the effective potentials Vu, Vu, V21 and V22 are defined by 

{p - P-){p - P+) 



V, 



11 



2 



16p5(p_ + po)(p+ + Po){p + Po)^(p^ + 2poP - P-P+ - (p- + P+)po) 
X [I6pi° + 128pop' + 32(p+(p_ + po) + po(p- + 15po))p' 

+8(f f9p|^ - 9(p_ + p^)pI - {IQpl + 25p+p_ + f6p2^)po - 16p_p+(p_ + p+))p' 

+ (1103p^ - 50f (p_ + p+)pI - (484p^ + 697p+p_ + 484p^)p2 

-196p_p+(p_ + p+)po + 288pV'+)p' 

+Po(764p^ - 959(p_ + p^)pI - (827p2_ + 634p+p_ + 827p^)p^ 

+325p_p+(p_ + p+)po + 1152pV+)p' 

+ (252p6 - 990(p_ + p+)pI + (-722p2_ + 193p+p_ - 722p^)p^ 

+ (p_ + p+)(72p2_ + 1159p+p_ + 72p^)p|; 

+3p_p+(f6p2_ + 611p+p_ + 16p^)p^ - 120p2_p^(p_ + p+)po - 96pip^)p" 

-2po(228(p_ + p+)pI + (53p2_ - 550p+p_ + 53p2^)p^ 

-(p_ + p+)(lllp'_ + 709p+p_ + lllp^)p;^ - 6p_p+(7p2_ + f04p+p_ + 7p^)p2 

+249p^_p^(p_ + p+)po + 180p-lp^)p-^ 

+ (-4(p_ + po)3(5po - 4p_)pt + Po(p- + Po)'(28p2_ - 559poP- + 151p^)p^ 

+pl{p. + po)(-12pi - 927pop'„ + 450p2p_ + 187p;^)p^ 

+P-Po(P- + Po)(-44p2_ - 2f3pop_ + 850p2)p+ 

V-Po(-20p'- + ISlpoP- + 187p2))p2 

-3po(p- + Po)(p+ + Po)(P-P+ + (p- + p+)po)(-16p^p+ - 7p_(p_ + p+)poP+ 

+3(3p2_ + 34p+p_ + 3p2^)p2)p 

+39p-p+Po(p- + Po)(p+ + Po)(p-P+ + (p- + P+)Po)^] , (51) 

13 



v^ 
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V, 



21 



2v^VP 


-P+{P- 


-P-){P- 


- P+) (P^ + 2poP - 


- P-P+ 


-(P- 


-+P+)Po) 








P^(p + Po)2 








v^VP- 


P+ip- 


-P-)ip- 


- P+) (P^ + 2poP - 


- P^P+ - 


-ip~ 


+ P+)po) 



2p4(p„+po)(p+ + Po)(p + Po)2 



(52) 



(53) 



V, 



22 



{P- P-)ip- P+) 



[I6p6 + 64poP^ + 96p2p4 



16p5(p_ + po)(p+ + Po)(p + PoY 

+8po(7po - (P- + P+)Po - P-P+)p^ 

+ (15p^ + ll(p_ + p+)pI + (12p2_ + 71p+p_ + 12p^)p2 

+60p_p+(p_ + p+)po + 48p2_P+)p' 

+5po(P- + Po)(p+ + Po)(16p-P+ + (p- + p+)po)p 

+39p-P+Po(P- + Po)(p+ + Po)] • 



(54) 



In order to discuss the stability for K = ±\ modes, we rewrite Eq. (H9|) and Eq. (l50|) in 
the following form 



- -^$ + V(p)$ = a;'$, 
dpi 

where $ and V represent the master variables and the potential given by 



(55) 



$ 






V 



Vn Vi2 
V21 V22 



(56) 



Multiplying both sides of Eg. (135]) by $^", and integrating it, we obtain 

2 



(ip* 



d^ 



dp^ 



+ $ty$ 



dp^ 



UJ 



dp* 1^1 , 



(57) 



I 1 9 + 

where |$| := #'#. Because the boundary term vanishes, if the second term of the integrand 
in the left hand side of Eg. fl57p is real and positive everywhere, there is no w^ < mode. 

We transform the term $"'" V$ as follows. We diagonalize the matrix V using the unitary 
transformation 



V^ll|^l|' + V^22|^2 



(58) 
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Here, V and ^ are defined by 

* = Q$=|""M, v = gvQt=|-^; |. (59) 





From Eg. ( 1571) and Eq.( l58|) . we can see that the sufficient condition for cj^ > is both Vu 
and V22 are real and positive in the domain — 00 < p* < 00. Here, Vu and V22 are written 
by using the components of matrix V as 

Vu + V22 + Viyn-V22y + WuV2, 
^11 - 2 ' 

Vu + V22 - ^{Vu-V22Y^W^2V2X ,^^. 

^22 - ^ • (bUj 

One can easily check that both Vu and V22 are real everywhere. Then, we split the condition 
both Vu > and V22 > to two conditions; (i) det V{= det V) > in the domain —00 < 
p=K < 00, and (ii) Vu > and V22 > at one point out side of the horizon. 

First, we show the condition (ii). Far from the horizon, Vu and V^2 behave as 

y^^ ^ 16(p- P-)(p-p+)p^° + C'('i')(6i) 

16p5(p_+po)(p^ + po)(p + po)3(p2^2pop-p_p+- (p_+p+)po)2 Vp/' 

V lQ{p-p-)ip-p+)p' (1\ . . 

'' i6p^{p^+po){p^+po){p+por v;- ^ ^ 

Since the dominant terms of (16T]) and (l62l) are positive, the condition (ii) is fulfilled at great 
distance. 

Next, we consider the condition (i). Typical profiles of det V are plotted in FigsJUlHl In 
the case a/3(5/(2M) < 1, because there are domains in which the det V^ is negative (see 
Figs|S]and[n]), we have not show the stability for this case. On the other hand, we can show 
the positivity of det V in maximally charged case \/3Q/{2M) = 1 (see FigH]). In maximally 
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charged case, we can write det V as 

256poVo(po - P+ + p)^(2po - P+ + p)' ^ ^ ^ 

+ (3072po + 1024p+)(p - p+)i2 ^ (i7408po2 + 12288p+po + 1536p^)(p - p+)" 

+ (59136po^ + 67840p+po2 + 18432p^po + 1024p:^)(p - p+)i° 

+ (132064pV + 219712p+po^ + 99552p^po2 + 12288p^po + 256p^)(p - p+f 

+ (203008po^ + 459968p+po^ + 309120p^po^ 

+65216p^Po2 + 3072ptpo)(p - P+)' 

+ (219328po^ + 656000p+po^ + 611232p^po^ 

+195648p:^Po^ + 1Q09QpXpo^){p - p+Y 

+ (165520po^ + 659024p+po^ + 819120p^po^ 

+366448p^Po' + 47104ptpo')(p - P+f 

+ (83521po^ + 472188p+po^ + 788454p^p~o^ 

+449980p^Po' + 84641p>V)(p - P+? 

+ (25572po^ + 234624p+po^ + 574056p^po^ 

+378400p;tPo^ + 96020ptpo^)(p - P+Y 

+ (3780po^° + 73616p+po'^ + 313560p^po^ 

+233104p:tpo' + 69188ptpo')(p - p+)' 

+ (11520p+po'° + 113280p>V + 110848p:tPo' + 32640pV~o^)(p - P+)' 

+ (19968p^Po'° + 36864p=^po' + 10752ptpo')(p - P+) 

+2048p~oVt + 6144po^Vl] , (63) 

where po is defined by Eq. (l57|) . From (1551) and fl^ . we can see the positivity of det l^ 
out side of the horizon. Since we have shown the condition (i) and (ii), we conclude that 
the background spacetime is stable against perturbation for K = ±1 modes in maximally 
charged case. 
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r„detV 




FIG. 4: The determinant of the potential matrix r^ det V in maximally charged case. 



400 



r:idetV 



3 Q/2M=2/3 




FIG. 5: The determinant of the potential matrix r^ det V in v3Q/2M = 2/3 case. 
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r„detV 




FIG. 6: The determinant of the potential matrix r^ det V in v3Q/2M = 1/3 case. 
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3. K = mode 



For the K = mode, we set hf^^ and 6Afj, as 



+2/i+_(p)e-'"V+a- + /i33(p)e-'"W3, (64) 

SAi.dx" = 6AA{p)e-''''dx^ + 6A3{p)e-''^'a\ (65) 

We set the gauge ^^dx^ as 

^f.dx'^ = Up)e-'^'dx^ + e3(p)e-*"V3, (gg) 

and under the gauge transformation, the metric perturbations transform as 

M p\p + Po) ' ^ ' 

htp -^ htp r- ^6 - lu^p + -Z-, (68) 

pip - P-){p - P+) dp 

, ^ , (p-+p+ + Po)p' + P-P+(-2p-po) rf^p .RQ^ 

p(p-P-)(p-P+)(p + Po) dp 

hfi -^ /its - «w6, (70) 

/i.- ^ /i.- + '^^"^-^^r'"y''^'°^ ep^ (72) 

P(P + Po) 

, ^ , , (P-P-)(P-P+)P0(P- +Po)(P++Po) ^ (^o\ 

pip + por 

So we can choose the condition ^ 

/i+_ = 0, /i« = 0, ht3 = 0. (74) 

In addition, for the gauge transformation ( 13T1) . the electromagnetic perturbations transform 
as 

Mi ^ Mi-ia;x(p), (75) 

M, ^M, + ^, (76) 

Ms ^ Ms. (77) 



Note that we can not choose this gauge condition for static perturbation. 
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We can choose the gauge condition 

5At = 0. (78) 

Substituting Eqs. (EH), (ESD, ([74]), and ([78]) into Eqs. ([21]) and ([28]), we obtain perturbation 
equations whose exphcit forms are given in Appendix [X] Fortunately, in these modes, the 
metric and the electromagnetic perturbations are decoupled. For the metric perturbations, 
the master equation becomes 

- J^^^OG + Vog{p)<^og = C^^^OG, (79) 

where the master variable $og(p) aiid the potential Vqg are defined by 

$0C(P) .- ^1/4(4^ + 3^^) ^33(P), (80) 

T. {P - P+){P - P+ + P-) 



16p5(po + p_-p^ + p)3(3(po + P-- P+) + 4p)^ 
X [(128po + 128(-p_ + 3p+))(p - p+f 

+ (528p^ + 32(-p_ + 55p+)po + 16(-35p^ + 98p+p_ + 17p^))(p - p+Y 

+ (720p|^ + 48(13p_ + 63p+)p2 + 16(-57p'_ + 330p+p_ + 79p^)po 

+ 16(-51p^ + 141p+p^ + 63p^p_ + 7p^))(p - p+)3 

+ (315p^ + 60(8p_ + Alp+)pl + (-450p^ + 6376p+p_ + 1850p^)pg 

+ (-1080pl + 5372p+p^ + 3056p^/5_ + 460p^)po 

-465pi + 35p^ + 328p_p^ + 1206p^p^ + 1456p^p+)(p - p+f 

+ ((792p+ - 81p_)p^ + 12(-27p^ + 239p+p_ + 98p^)p|^ 

+ (-486pi + 3852p+p^ + 2834p^p_ + 520p^)p^ 

+4(-81pi + 567p+pi + 535p^p^ + 209p^p„ + 18pt)po 

+p_(-81pi + 492p+pi + 482p^p^ + 316p^p_ + 71pt))(p - p+) 

+36p^p+(3p_ + 8p+) + 12pj^p+(36p^ + 73p+p_ + 16p^) 

+4p2p+(162pi + 225p+p^ + 121p^p_ + 8p=^) 

+4p^p+(27p;L + 3p+p^ + 25p^p_ + 9p^) 

+4pop_p+(108pi + 81p+p2_ + 98p^p_ + 17p^)] . (81) 
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The master equation for the electromagnetic perturbations is given by 

dpi 
where ^oe{p) and Vqe are defined by 



%e{p):=p'^\p + Po?^'SA,{p), 



(82) 



(83) 



Vn 



OE 



[(8po-8p- + 24p+)(p-pH 



(P-P+)(P-P+ + P-) 
16p^(po + P- - p+ + p)^ 
+ {15pl + (2p_ + 50p+)po - I3pl + Q3pl - 2p_p+)(p - p+)2 

+ (5(-p_ + 8p+)pg ^ ^Q(_ ~2_ _ 3^^ ~_ ^ ^2p^)po 

+p_(-5p2_ - 70p+p_ + 123p^))(p - p+) 

+4p2_p+(15p+ - llp_) + 4pop_p+(31p+ - 22p_) + p2(64p^ - 44p_p^ 



)]• (84) 



From fl38|) . f lHT]) and fl84|) . we can see Vqg > and Vqe > out side of the horizon. The 
stabihty has been shown against the perturbation for K = mode. Typical profiles of the 
potentials are plotted in FigsJTllTOl 
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FIG. 7: The effective potential Vqg in maximally charged case. 



B. K = =b( J + 2) modes perturbation 

As noted in the previous section, the highest modes for /i++ and h are always decoupled 

for arbitrary J. Since these are gauge invariant, it is straightforward to get the perturbation 
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FIG. 8: The effective potential Vqg in V^Q/2M = 1/2 case. 
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FIG. 9: The effective potential Vqe in maximally charged case. 
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FIG. 10: The effective potential Vqe in V^Q/2M = 1/2 case. 
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equation for h^^^ whose explicit form is given in Appendix |X1 Defining the new variable 
we obtain the master equation 

- j-2^j + yAp)^j = ^^^j^ (86) 

where the potential Vj{p) is defined by 

T/ {P - P+){P - P+ + P-) \^R~iln I j\2J2 

lQpo{po + p-)p^{po + P-+p-p+r'- 
+4plp43p- + (57 + 60J + 16 j2)p+)po + 4pV+(9p- + 2(37 + 42J + 12 j2)p_^)p2 

+4p_p+(9p_ + (41 + 52J + 16J^)p+)pI + 4p+(3p_ + 4(2 + 3 J + J')p+)p^ 

+ (32p3_(2 + J)V+(P- + 2p+) + p^(-9p^ + 2p_(249 + 240J + 64j2)p+ 

+ (711 + 736J + 192J^)pI)po + p„(-27p^ + 4p_(175 + 168J + 48j2)p+ 

+ (655 + 704 J + 192 J2)p^)p2 + (-27p^ + 2p_(209 + 208 J + 64j2)p+ 

+8(25 + 28J + 8J')p^)p|^ + (-9p_ + 8(11 + 12J + 4j2)p+)/5^)(p - p+) 

+ (16p^(2 + jy{pl + 8p_p+ + 6p^) + p-(p^(207 + 240J + 647^) 

+2p_(735 + 736J + 192j2)p+ + (739 + 752J + 192j2)p^)po 

+ (p^(257 + 336J + 96J2) + 8p_(175 + 176J + 48J^)p+ 

+ (355 + 368 J + 96 J^)p^)pg + (p_(149 + 208 J + 64 J^) 

+2(221 + 224J + 64j2)p+)p|; + (35 + 48J + 16J')p^)(p - p+)' 

+ (64p_(2 + jy{pl + 3p_p+ + pD + 8(p^(87 + 92J + 247") 

+p_(187 + 188J + 48j2)p+ + 8(2 + J)V+)po + 8(8p_(10 + IIJ + 3J^) 

+ (91 + 92J + 2AJ^)p+)pI + 8(25 + 28J + 8j2)p|^)(p - p+f 

+ (16(2 + J)'(6p^ + 8p_p+ + pI) + 16(2 + J)(p-(23 + 12J) + 8(2 + J)p+)po 

+16(22 + 23J + QJ^)pl){p - p+Y 

+ (32(2 + J)2(2p_ + p+) + 64(2 + Jfpo){p - P+f + 16(2 + J)\p - p+f] . (87) 



Since we can see Vj > out side of the horizon from fl38|) and flHTI) . the stability has been 
shown against perturbation for K = ±( J + 2) modes. 
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IV. SUMMARY AND DISCUSSION 



In this paper, we have extended the stabihty analysis for perturbations of the SqKK 



black holes given in 56| to the SqKK black holes with charge 3l|. We have expanded the 
perturbation variables in terms of Wigner function -D^jv/ which have three indices J, M, K, 
and derived the master equations for ( J = 0, M = 0, K = 0, ±1, ±2) modes which have 
SU{2) symmetry and the highest modes {K = ±(J + 2)). Using the effective potential 
functions in the master equations, we have discussed the stability for these modes. 

Since the modes oi J = 0,M = 0,K = ±2 and the highest modes contain only the 
gravitational perturbation, the perturbation equation for each mode can be reduced to a 
single master equation. In the mode of J = 0, M = 0, -ft' = 0, the perturbed equations 
can be reduced to two decoupled master equations. As for the modes of J = 0, M = 
0,K = ±1, we have obtained the coupled master equations for the gravitational and the 
electromagnetic perturbations, in contrast to the case for Reissner-Nordstrom black holes of 
which these perturbations are decoupled against all modes [l8| . We guess that the coupling 
for the gravitational and the electromagnetic perturbations comes from the differences of 
the asymptotic structures. 

We have shown from the positivity of the effective potential functions of the master 
equations that the modes of J = 0, M = 0, -ft' = 0, ±2 and K = ±( J+2) are stable as those in 
the neutral case. In the case ofJ = 0,M = 0,-ft' = ±l, the effective potential functions form 
a matrix. In maximally charged case, we have shown the stability for J = 0,-ft' = ±lby using 
the matrix of effective potential functions. In general, because there are domains in which 
the effective potential functions are neg ative, we must find an appropriate S-deformation 



function similar to the discussions in 
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18| to show the stability for J = 0, M = 0, iT = ±1 
modes. Even in neutral cases where the gravitational and the electromagnetic perturbations 
are decoupled, the S-deformation function have not been found yet. It is hard to show the 
stability of these modes in the case of charged black hole. To show stability for J = 0,K = 
±1 modes in general is open issue. 



Since the stability for J = modes is suggested in the limit of both neutral case |56l.l57| 
and extremal case, and the instability empirically appears in the lower modes, we expect 
SqKK black holes with charge are stable. 
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Appendix A: the perturbation equations 

In this section, the precise form of the perturbation equations are given for each mode. 

1. K = ±2 modes 

In K = 2 mode, there is only (++) component of linearized Einstein equation which is 
given by 



■[4p^ + 16pVo - 2pop+p_(po + P+)(po + p- 



h++ 



V(P + Po)^(po + P+)(po + P-) 

+V(5po - P+P- - Po(p+ + P-)) + PPo(Po + P+)(Po + P-)(-6p+P- + Po(p+ + P-)) 

+V(3po + P+P-{p+ + P-) + Po(p+ + P+P- +P-))+ p'(2p^ - 6p>'_ + pl{p+ + p_: 

-3poP+P-(p+ + P-) + Po(3p+ + 3p^ - 2p+p_))] 
2p^ + 2pop+p_ - 3p^(p+ + p_) - p(-4p+p_ + po(p+ + p_)) dh++ 

2p2(p + po)2 dp 

{p- P+){p- P~)d^K+ P^ 



., rW /l++ = 0. (Al) 

2p(p + po) dp^ 2(p-p+)(p-p_) ++ ^ ' 

The equation for K = —2 mode can be obtained by just replacing /i++ to h in Eq. (lAip . 

2. K = ±1 modes 

In K = 1 mode, there are (t+), (p+) and (+3) components of linearized Einstein equation 
and there is (+) component of Maxwell equation for perturbations which are given by 



SGt+ - 26Tt+ 



[p5 + 3pop^ + 3py + (pI - 2(p_ + p+)pI - 2(p_ + 



2p3(p_ + po)(p+ + Po)(p + Po) 

P+ypo - 2p_p+(p_ + p+))p'^ - (p_ + Po)(p+ + Po)((P- + P+)p 
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-^P~P+)P + P-P+PoiP- + Po)(p+ + Po)]ht^ 



(p- P-){p- P+)podht+ 
2p2(p + po)2 dp 



{p - P-){p - P+) d'^ht+ iu{p- p_){p- p+) iuj{p - p_){p - p+) dhp+ 



2p{p + po) dp^ p^ip + po) ''^ 2p(p + po) dp 

V3^p-p+{p- - p){p+ - p) d5A^ 



P^ip + Po) dp 



0, (A2) 



SGp+ - 25Tp+ 

(p^ + 2pop - p_p+ - (p_ + p+)po)^ ^ ia;p(2p + po) ^ 

2p(p_ + Po) (p+ + po){p + PqY ^^ 2(p - p_) (p - p+) (p + Po) 

iuV^^/pZp^ iup^ dht+ u^p^ ^ 

:^^+ + 7^7 ^7 ^—^ TT7 77 7 V = 0' (A3) 



{p-p-){p-p+) '^{p- P-){p- P+) dp 2(p-p_)(p-pH 



5G+3 - 2(5T+3 

■ [(2p5 + (-p_ - p+ + 5po)p' + 2(3p2 + p.p+)p^ - (5(p_ + p+)pI 



2p2(p + po)4' 
+ (3p^ + 8p+p_ + 3p^)po + 3p_p+(p_ + p+))p^ + 2p_p+(p2 + 2(p_ + p+)po + 2p_p+)p 



^P (P + 2poP - P~P+ - (P- + P+)pi 



•o) 



+P-P+Poip-P+ + (P- + P+)po))]hp+ + -^ ^7 r-, ■ r^ ht+ 

oi^ - ^'- - ){p + poy 

0, (A4) 



2(p-p_)(p-p+)(p + po)2 
Kp - P-){p - P+){P^ + 2poP - P-P+ - (p- + P+)po) dhp+ 



2p(p + po)3 dp 

-1 6A+ ((p_ + P+ +po)p2-2p_p+p-p_p+po)rfM- 



sm 6* p2(p_ +po)(p+ +Po) p^{p + po)^ dp 

UJ^P . . V^^/P^P+{p- - P) (P+ - P) (2p + Po) , 

-oA+ ^^77 TT^ ht+ 



(P - P-)(P - P+)(P + Po) 2p3(p _ p_)[p - p^)[p + po)3 

V3VP^P+(P- - P)(P+ - P) dht+ iujV3,/p7p^{p^ - p) (p+ - p) 
+ 2p2(p-p_)(p-p+)(p + Po)2 dp +2p2(p_p„)(p_p^)(p + po)2 ''+ 

_ (PZ_Pz)(qP±) f^] = 0. (A5) 

p2(p + po)2 dp^ ' 

The equations for if = — 1 mode can be obtained by taking Eq. (]A2p . (lA3p . (]A4p . and 
( lASp to its complex conjugate equations, and using the relations Jia^ = hA-, h^s = h-s and 
6A. = 6A_. 
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3. K = mode 

In i^ = mode, there are (tt), (tp), (pp), (t3), (p3), (H — ) and (33) components of lin- 
earized Einstein equation and there are (t), (p) and (3) components of Maxwell equation for 
perturbations which are given by 

SGtt - 2STtt 
- '-'-^' - '-^^' - '-^ h, + ^^^^47^^^^^^^[4p^ + 4(p_ + p, + 3po)p^ 



4p6(P-+Po)(p+ + Po) 4p6(p + po) 

+ (6po + 2(p- + p+)po - 9p_p+)p2 - 10p_p+poP - 3p-P+pl]hpp 
(P - P-)(P - P+)(4p-'^ + (Po - 3(p^ + p+))p^ + Ip^P+P - P-P+Po) dhs3 
4p5(p_+po)(p+ + po)(p + Po) dp 

(p-p_)3(p-p+)3(4p + 3po)rf/ipp {p-p^y{p-p+y d^hs 



+ 

+ ^^^^^^'T.^::r^^^'^^ M, = 0, (A6) 



4p5(p + po)3 dp 2p4(p_+po)(p+ + Po) dp^ 

ibJ^n^fp Zpl{p^ - p) {p^ - p) 



6Gtp - 26Ttp 

-iujp - p+){p - p-)(3po + 4p) ia;(po(p+p_ - p^) - p{p_p + p+{p - 2p_))) 

4p2(po + p)2 '' 4(po + p+)(po + P-)(p+-p)(p--p)p2 '' 

^^uiPo_U^dh^^^^ 

2(Po + P+)(Po + P-)p dp 



SGpp - 26Tpp 

{p- - p)po{4:p + Po) + P+(Po(po - 2p-) + p{p^ + 4po)) 



4p2(p - p_)(p - p+){p^ + po)(p+ + Po)(p + Po) 
-4p2 - 8pop - 3pg + P-P+ + (p- + p+)po 

4p(p + po)3 
4p2 + (2po - 3(p_ + p+))p + 2p_p+ - (p_ + p+)po c?/i33 
4p(p-P-)(p-p+)(P-+Po)(p++Po) dp 

iupjAp + 3po) ^V(p + Po)'^ , 

'2(p-p-)(p-p+)(p + Po) *' 2(p-p_)2(p-p+)2(p_+po)(p^ + po) ■'' 
iuVSy/pZp^ 



(p - P-)ip - P+) 



6Ap = 0, (A8) 



SGt3 - 25Tt3 

V^yp^P+JP- - p){p+ - P) d5A3 iujp - p^){p- p+) 
P^{p + Po) dp p2(p + po) 
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hp3 



iuj{p-p_){p- p+)dhp3 



2p(p + Po) dp 


) 


SGp3 - 2STp3 




iu^/li^p^p+ ^ ^ 


wV' 


{p-p^){p-p^) ' ■' 2(p- 


-P-)ip-P+) 



(A9) 



h 



p3 



0, 



(AlO) 



5G+- - 26T+^ 

2poP^ - (Po + 5(p- + p+)po + P-P+)p^ - 2po{{p- + p+)po - 3p_p+)p + 4p_p+p2 



1 



2p2(p + Po)3' 



2p3(p_ + po)(p+ + Po)(p + Po) 
-2(p_ + p+ + po)p" + (p'_ + 5p+p_ 



h. 



33 



pI 



3p^-2(p_+p+)po)p' 



+ (2(p_ + p+)p2 + 2(p2„ + 5p+p_ + p2^)po - P-P+{P- + P+))P' 
+P-P+(Po - 4(p- + P+)po - P-P+)p - P-P+ip- + P+)pl]hpp 



P' + ((P- + P+)Po - P-P^)P - 2P-P+P0 fc ^ 2.c.y3VpIpTM, 



P^(P-+Po)(p+ + Po) 
(P - P-)(P - P+)(-P- - P+ + 2p) rf/ipp 



(ip 



2p(p + Po) 
wV(P + Po)' 



+ 



{p - P-){P - P+){P + Po) d'^h: 
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P(P- +Po)(p+ + Po) dp'^ 



{p- P-)ip- P+){p- + Po)(p+ + Po) 



dp 



dp 



iojp{2p^ + (-p_ - p+ + 3po)p^ - 2(p_ + p+)pop + p-p+po) 
(p-P-)(p-P+)(p + Po) 



tp 



(All) 



SG33 - 25T33 

- ^''"V^p + Pol^''"^ ^^^' + 12poP^ + i-5pl - I3p^p- - 3pl + 3pl - 12(p_ + p^)po)p 
+ (-2(p- + P+)Po + 2(p'_ + 3p+p_ + p'+)po + llp-p+(p- + p+))p^ - P-P+U 



-2(p- + p+)po + 7p_p+)p + p_p+po((p- + p+)po - 2p_p+)]hpp 
(p-P-)(p-P+)(P-+Po)(P+ + Po)r/ ,2 



4p(p + Po)^ 

,dhp 
dp 
iup{p_ + po){p+ + Po) 



o < { , \ M-'VP , 3(p_ + po)(p+ + Po) 

-2p_p+ + (p_ + p+)po)]-^ + - 



[(-4p^ + (3(p_+p+)-2po)p 
^33 



2(p-P-)(p-p+)(p + Po)=' 



4p(p + po)3 
[(4p=' + 3(-p_ - p+ + po)p' - 2((p_ + p+)po - p-p+)p 



iuV3^p^p+{p- + po)(p+ + Po) . . 
+p_p+Po)J V + -— — -^ 6Ap 

i(^p'^ip-+ po)ip++ Po)dhtp a;V(P- +Po)(p+ + Po) 



(P + Po)^ 



dp 



2(p + Po)2 



pp 



0, 
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(A12) 



-1 
V(P + Po)-'(po + P+)(po + P-) 
[^VP-P+PoiP + Pofhs + V3y/p-p+{p- + po){p+ + Po)((P- +P++ Po)p^ 

-2p_p+p - P^P+Po)hpp - V3VP-P+p(p + /^o)^^^ + ^VP-P+(P- " P)(P+ " P)P 

(P- + Po)(p+ + Po)(p + Po)^ + t^Sp\p- + Po)(p+ + Po)(p + Pof^Ap 

+iujAp\p_ + po)(p+ + Po)(p + Po)'^] = 0, (A13) 



-1 



V(P + Po)'(Po + P+)(Po + P-) 

[~iuj^/?>^p^p+{p + po)^/i33 + icjV3^p_p+(p_ - p)(p+ - p)(p_ + po)(p+ + po)hpp 

-w24p3(^_ ^ p^^^p^ ^ p^^^p ^ p^)SAp] = 0, (A14) 

iuVS^^pZp^ 1 

2p2(p-+Po)(p+ + Po) ^' p2(p + Po)2 ' 
^ V + (-P- - P+ + Po)p^ - P-P+Po dSAs ^ {p-p_){p-p+) d^6As 
p3(p_+po)(p++po)(p + po) dp p2(p_+po)(p+ + po) rfp2 

^'^(^ + ^°) M3 = 0. (A15) 



(P - P-)ip - P+)(P- + Po)(p+ + Po) 
4. K = ±(J + 2) modes 



In K = (J + 2) mode, there is only (++) component of linearized Einstein equation which 
is given by 






■[(J + 2) V + 4(J + 2)Vop' + P^((2J + 5)(3J + 4)p2 



2p3(p_+po)(p+ + po)(p + po)3' 

-(J + 4)(p_ + p+)po - (J + 4)p„p+) + 2p3((J + 2)(2J + 3)pI - J{p^ + p+)pl 

+ (2p2_ - (J - 2)p+p_ + 2p^)po + 2p_p+(p_ + p+)) + p\{J + 1)(J + 2)p^ 

28 



-(J - l)(p_ + p+)pI + (3p2_ - (J + 2)p+p_ + -ipDpl - -ip-p+ip- + p+)po - Qplpl 

+PPo{p- + Po)(p+ + Po)((P- + P+)Po - 6p_p+) - 2p_p+p^(p_ + po)(p+ + Po)] 
2p^ - 3(p, + p+)p'^ - ((p- + p+)po - 4p_p+)p + 2p_p+po dh++ 

2p2(p + po)2 rfp 

{p- P-){P- P+)d^K+ P^ 



,, -u;^/i++ = 0. (A16) 

2p(p + po) dp^ 2(p-p_)(p-p+) ++ ^ ^ 

The equation for K = —{J + 2) mode can be obtained by replacing h^^ to h in Eq. 

dAH]). 
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